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ABSTRACT 

An isolated, initially cold and ellipsoidal cloud of self-gravitating particles represents a 
relatively simple system to study the effects of the deviations from spherical symme¬ 
try in the mechanism of violent relaxation. Initial deviations from spherical symmetry 
are shown to play a dynamical role that is equivalent to that of density fluctuations in 
the case of an initially spherical cloud. Indeed, these deviations control the amount of 
particles energy change and thus determine the properties of the final energy distribu¬ 
tion, particularly the appearance of two species of particles: bound and free. Ejection 
of mass and energy from the system together with the formation of a density pro¬ 
file decaying as p{r) ~ r -4 and a Keplerian radial velocity dispersion profile, are the 
prominent features similar to those observed after the violent relaxation of spherical 
clouds. In addition, we find that ejected particles are characterized by highly non- 
spherical shapes, whose features can be traced in the initial deviations from spherical 
symmetry that are amplified during the dynamical evolution: particles can indeed form 
anisotropic configurations, like bars and/or disks, even though the initial cloud was 
very close to spherical. 

Key words: galaxies: formation; galaxies: elliptical and lenticular, cD; methods: 
numerical 


1 INTRODUCTION 

The collapse of an isolated, cold and initially spherical cloud 


the literature 

Henor 

1973l:IVan Albadal 1982l;lAarseth et all 

19881; Bertin 


200C 

: 

Boilv et all 20021; 

Jovce et al. 2009bl; 

Svlos Labinil 

2012j, 

20131) as it is considered a paradigmatic 


example of the violent relaxation mechanism. This is the 
main physical process through which a self-gravitating sys¬ 
tem reaches a quasi stationary virialized state far from ther¬ 
modynamic equilibrium. Numerical experiments have shown 
that the global collapse of the system occurs in a relatively 
short time scale r c ~ 1/y/Gpo, where po is the cloud mass 
density. For t < r c the violently changing gravitational field 
of the system considerably varie s individual particle energy 
dHenonlll964 iLynden-Bellll 19671 ) . This mechanism of energy 
exchange was called violent relaxation in order to distinguish 
it from the collisional t wo-body relaxation which occurs on 
much longer time scale (Binne v fe Tremai ne 201 ll ). 

One of the most peculiar features occurring in the 
collapse of a spherical and cold cloud is represented by 
the ejection of a relevant part of its mass (i .e., ss 30%) 


and e n ergy. While seve r al au t hors (see, e.g. VmiAlbada| 


( 1982jjj_ David fe Theunsl (1198911 : Aguilar fe MerrittF 1990jb 


iTheuns fe Davidl 1 19901) : iRov fe Pere d ( 20041) : Barnes et al.l 
( 20091) ) noticed this ejection, the importance of this process 


with respect to the mechanism of violent relaxation and to 
the formation of the virialized structure h as been overlooked . 
In a series of works djovee et al.ll2009bl lal; ISvlos Labinill20f^ . 
120131) we have shown that the ejection of mass and energy 
is indeed the crucial mechanism for the formation of the 
virialized structure density and velocity profiles. 

Mass and energy ejection depends on both the ini¬ 
tial shape of the dens i ty pr ofile and the initial virial ra¬ 
tio dSvlos Labinil l20l3 . I 2 OI 3 I ). Namely, the largest system 
contraction occurs for a cold and uniform initial density dis¬ 
tribution; if the density profile is described by a power-law 
function of the distance from the origin, i.e. p(r) ~ r _ “, 
the collapse gets softer when the exponent gets steeper. 
Eventually, for a > 2 the system does not contract vi¬ 
olently anymore and its final confi guration does not sub¬ 
stantially differ from the initial one dSvlos Labinill2013l ). On 
the other hand, when the initial velocity dispersion is large 
enough, i.e. b = 2K/W with b < —1/2 (where K/W is 
the initial kinetic/potential energy), random motions of the 
particles oppose to the system gravitational contraction, so 
that the cloud only slightly changes its initial configuration 
nn os Labin3l2012l ). 

Only when particles ejection occurs, the density and 
the radial velocity dispersion profiles of the remaining viri¬ 
alized structure show the same behavior: in particular it is 
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observed that the latter scales as p(r) ~ r -4 while the for¬ 
mer scales as ( v (!) ~ r _1 . If ejection does not take place then 
the power-law tails do not form. This tail is indeed made of 
bound particles with energy close to zero, i.e. with velocity 
only slightly smaller than the escape one. The universality 
of the density profile thus originates from the large spread 
in the energy distribution formed during the violent collapse 
of the cloud. Because of such large energy change, some of 
the particles can gain enough kinetic energy to escape from 
the system while other remain bound to it but with energy 
very close to zero. 

Therefore a spherical and cold cloud represents a very 
simple initial configuration showing a rich variety of non 
trivial dynamical behaviors during its evolution. This is how¬ 
ever a too idealized situation for an interesting physical ini¬ 
tial condition. As a further step towards a more realistic sit¬ 
uation we consider in this work the collapse of a cold and iso¬ 
lated ellipsoidal cloud of self-gravitating particles. We show 
that the initial deviations from spherical symmetry, finely 
tuned in numerical experiments, can be interpreted to play 
the same role of density fluctuations for the collapse of a 
spherical system, determining both the ejection of a frac¬ 
tion of the system mass and the formation of the density 
and radial velocity profiles decays respectively as p(r ) ~ r ~ 4 
and (Vr) ~ r^ 1 . In addition, we show that both the virialized 
structure and the ejected particles are characterized by non- 
spherical shapes, whose characteristics reflect, in a different 
way, the imprint of the deviations from spherical symmetry 
characterizing the initial clouds. 

The paper is organized as follows: we first describe in 
Sect[2]the initial conditions we used and the details of the 
simulations we have performed.Then in Sect[3] we discuss 
the main features of the violent relaxation mechanism for 
an initially ellipsoidal cloud. Post-collapse virialized parti¬ 
cles shows the same density and radial velocity profiles that 
are independent of the initial cloud shape: these features, to¬ 
gether with the strongly anisotropic configuration of ejected 
particles, are discussed in Sect|4] Finally we draw our main 
conclusions in Sect [5] 


2 INITIAL CONDITIONS AND SIMULATIONS 

We first describe the properties of the initial conditions 
that we considered for our controlled numerical experiments. 
Then we discuss the technicalities of the simulations. 


2.1 Initial conditions 

Initial conditions are represented by an isolated cloud of N 
particles of equal mass m. The total mass of the system 
M — m ■ N is kept constant and thus when the number 
of particles N is changed their mass scales accordingly, i.e. 
m = M/N. Note that the thermodynamic limit corresponds 
to 


lim m ■ N = M = const. (1) 

N—>oo-,m—>0 

Particles are distributed randomly with uniform density, i.e. 
density fluctuations are Poissonian and decay as AN/N ~ 
N~ 1/2 . 


Name 

c 

T 

<t> 

Sphere 

0 

- 

0 

Prolate 

X 

1 

ss x/2 

Oblate 

X 

0 

RS x/2 

Triaxial 

2x/(2 + x) 

1/2 

2x/(4 + 3x) 

Disk 

1 

0 

1/2 

Tiny Cylinder 

> 1 

1 

« 1 


Table 1. Oblate (a2 = (13 = ai -\-x > a±), prolate ( 0,3 = 0,2 + x = 
a\+x ) and triaxial ( 0,2 = (ai+as)/2, where as = a\-\-x) ellipsoid: 
in all cases we take 0 < x 1. 


The cloud has initially an 


ellipsoidal shape 0 



( 2 ) 


where we take hereafter 03 ^ 02 ^ an. The three eigenvalues 
of the inertia tensor can be calculated as 


Ai = ^M(a] + a 2 k ) (3) 

where i ^ j ^ k and i,j,k,= 1,.., 3: from the definition 
of the semi-axes we have Ai > A 2 ^ A 3 — note that Ai is 
oriented in the direction of the smallest semi-axis ai and A 3 
in the direction of the largest one 03 . 

To characterize the structure shape we define three dif¬ 
ferent linear combinations of Athe flatness parametei 0 


^ 1 , 

L ~ A 3 

(4) 

the triaxiality index 


A2 — A3 

(5) 

T = - 

Ai — A3 

and the disk parameter 


4> = Al ” As . 

A2 + A3 

( 6 ) 


Note that the combination of 1 , r and </> allows one to distin¬ 
guish not only between different type of ellipsoids (prolate, 
oblate and triaxial) but also between other shapes like bars 
and disks. In Tab[l]we report their values for some paradig¬ 
matic cases. 


2.2 Simulations 

We have performed N r = 10 different simulations of the 
same cloud but with different realizations of the initial Pois¬ 
son noise. We have computed averages of the various phys¬ 
ical quantities over the different realizations of the cloud. In 


1 Our units are such that, for the spherical case aj = a 2 = 03 = 
R c we have po = Mj (4 tt R/. /3) = 1 gr/cm 3 so that the character¬ 
istic time scale for the collapse (see Eq[ 8 j) is r c = 2100 seconds 

2 Note that these parameters are generally defined in function of 
the semi-axis ai, a 2,03 and not in function of the eigenvalues (see 
EofTli. However for small deformations of a perfect sphere (those 
that we considered here) these definitions are almost equivalent. 
The advantage of using the eigenvalues to define these parameters 
is that this choice allows one to describe more easily shapes that 
differ from ellipsoids. 
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addition, we considered two set of clouds with N = 10 4 and 
N = 10 5 particles and same total mass. 

We have used the parallel version of the publicly avail - 
able tree-code Gadget dSpringel et al.ll 200 ll : ISpringelll2005l l 
to run N-body simulations. We used a softening length e, i.e. 
the scale below which the two-body potential is not New¬ 
tonian anymore, such that e/£{t = 0 ) -C 1 , where £(t = 0 ) 
is the initial inter-particle distance. In order to choose a 
proper value for e that does not affect the gravitational dy¬ 
namics during the whole time range of the simulations (i.e. 
t < 10 t c ), we have considered a necessary (but a priori not 
sufficient) condition: e must be smaller than the system size 
at all times. As typical radius of the system we took the 
gravitational one: 


Rg(t) 


GMb(t) 

\W b (t)\ 


(7) 


where M b (t) and 114(f) are respectively the mass of the 
bound system and its potential energy at time t. We numer¬ 
ically determine the minimal value R™ ln of R g (t) during 
the time evolution, that occurs for times close to the col¬ 
lapse time T c ~, we thus we require that R™ ln e. Clearly at 
t « t c the inter-particle distance is £{t c ) -C £, so that two- 
body scatterings can highly differ from a pure Newtonian 
situation. However the typical time scale of two-body scat¬ 
terings is far longer than the typical collapse time scale: for 
this reason collisional effects should not be dominant during 
the collapse phase. 

In order to test the effects of collisionality we have per¬ 
formed several tests by varying the smoothing length e in 
Gadget (and the appropriate time-step parameters) and we 
have carefully foll owed the co l lapse. Our conclusion, that 
agrees with that of I Joyce et alj d2009b|j . who considered the 
collapse of a initially cold, spherical and uniform cloud (for 
which the collapse is stronger), is that as long as we have 
e < Rmin, collisionality does not play an important role in 
the violent collapse phase of these systems. 

In addition, we found that the necessary condition 
e <C R™ ln is also a sufficient one, as not only the numer¬ 
ical integration conserves total energy and momentum up 
to a few percent but also macroscopic quantities describing 
the system we are interested in, as the density and veloc¬ 
ity profiles, the global shape, etc., do not show a detectable 
dependence on e. In particular, we have found a good con¬ 
vergence of the results for R™ ln /e € [10, 200]. On the other 
hand when £ ^ R™ ln the gravitational dynamics is modified 
at the scale of the system and one detects large deviations 
with respec t to the £ < R™ m case (see also discussion in 
I Joyce et alj d2009bl l : |Svlos Labinil d2012l Yh 

Beyond the softening length, the precision of a Gadget 
simulation is determined by the internal time-step accuracy 
r/ and by the cell-opening accuracy parameter of the force 
calculation of- We have chosen the time-step criterion 0 of 
Gadget with p = 0.025. In the force calculation we employed 
the (new) Gadg et cell opening criter i on with a force accuracy 
of cxf = 0.001 (iSpringel et aljfiooll : ISpringelll2005l '). 



Figure 1. Example of the behavior of the gravitational radius 
(scaled to its initial value R g ( 0)) as a function of time when the 
initial condition is a cold spherical cloud. It is also shown the 
method we choose to estimate the duration of the collapse. 


3 VIOLENT RELAXATION MECHANISM FOR 
AN ELLIPSOIDAL CLOUD 

We first review in Sect l3.li the main elements of the violent 
relaxation mechanism for initially spherical clouds with dif¬ 
ferent discretization but same total mass. Then in Sect l3.2l 
we discuss the case of an initially ellipsoidal cloud consider¬ 
ing how to map one problem into the other. 


3.1 The spherical case 

For an idealized initial spherical cloud all particles have the 
same collapse time, i.e. they arrive at the center simultane¬ 
ously: 


37T 


32Gp 0 


(8) 


where po = 3M/(4nRc). In the thermodynamic limit (i.e., 
when fluctuations can be neglected) the collapse can be de¬ 
scribed as the simple gravitational contraction of a perfectly 
homogeneous cloud. This simple idealization well describes 
the observed collapse for t < r c . Because of the effect of den¬ 
sity fluctuations, the same a pproximation fails t o follow the 
system behavior for t ^ r c (I Joyce et al.l I2009bl '). i.e. when 
the system size should nominally reduce to zero to then re¬ 
expand in a periodic way. Instead, for t > r c a simulated 
cloud relaxes into a quasi stationary state with its size re¬ 
maining almost constant (see FiglT]). 

_ Several auth ors (see, e.g ., lAarseth et al.l (1 1988ll : 

iBoilv et~ahl (l2002fl ; I Joyce et al.1 d2009bl ') and references 
therein) found that the minimal gravitational radius of the 
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N 

Figure 2. Gravitational radius (see EqO for an initial spherical 
cloud with different values of N. 



Figure 3. Behavior of the spread of free fall times AT (normal¬ 
ized to r c ) for an initial spherical cloud and different values of 
N. 


bound structure R™ in scales as 


R 


min 

9 


N~ 1/3 . 


(9) 


This relation physically means that the system stops collaps¬ 
ing, reaching its minimal size, when N density fluctuations 
on the scale of the systems itself go nonlinear: Eq[9] nicely 
fits the observed behavior as shown in Fig(2] 

When N is finite, the free fall time spread is related to 
the initial density fluctuations. From Eq[8]we simply get 


AT 

Tc 


^ ~ AT 1 ' 2 , 

po 


( 10 ) 


which approximates well the data (see Fig[3]). Note that in 
the simulations AT has been estimated as the width of R g (t) 
around its minimum @. 

During the collapse a fraction of the particles may gain 
enough kinetic energy to escape from the system. In par¬ 
ticular the fraction of ejected particles shows a slow log- 
depende nce with N (see Fi g l4)l in agreement with the re¬ 
sults of I Joyce et alj ( 2009bl) . We recall that there was no 
theoretical explanation for this behavior and we will come 
back on this point in Sect l3.2l The mechanism of particles 
energy gain originates from the coupling of the growth of 
density fluctuations with the system finite size. Particles 
originally placed close to the system boundaries develop a 
net lag with respect to the bulk, thus arriving at the sys¬ 
tem center when the others are already re-expanding. These 
particles, by moving in a rapidly varying gravitational field 


3 To estimate AT we define the two values of the gravita¬ 
tional radius for t < t c and t > r c such that R* = + 

( Rg Syn — R™ ln )/2 where Rg ayn is gravitational radius reached 
for t > T c - Then the duration of the collapse is defined as 
A T = t 2 (R*)-t 1 (R* g ) (seeFigO. 


can gain kinetic ener gy and can drastica lly change the whole 
energy distribution (I Joyce et alj|2009bl ). 

The energy change particles undergo when passing 
through the center has also an angular dependence that re¬ 
flects the initial deviation from perfect spherical symmetry. 
Indeed, the flatness parameter (Eq(4j is initially different 
from zero because of Poisson fluctuations 

i(0) ~ ' (11) 

This initial small deviation from perfect spherical symmetry 
is amplified during the collapse phase because of the large 
particles energy change. Eg II 1 1 implies that the system is 
not initially perfectly spherical, but it can be described as a 
ellipsoid with Ai > A 3 (i.e. ai < 03 ). According to the en¬ 
ergy gain mechanism we described above, particles arriving 
later at the center, i.e. those initially placed in the range 
d 2 < r < 03 , get the largest energy kick. Some of them 
are ejected but others are still bounded for t > r c . A sim- 
ple ansatz to desc r ibe th e deformation was introduced by 
Worrakitpoonponl (|2015i l: ISvlos Labini. Benhaiem fe Jovcel 
( 20l4) : if AT is the time during which energy is exchanged, 
we can estimate that the difference between the major and 
minor semi-axis is 


b{t > t c ) ~ AT ■ Av c , ( 12 ) 

where v c is the typical particle velocity at the collapse time. 
This simply means that difference in the linear dimension 
between 03 and ai < 03 of the structure depends on the time 
interval AT characterizing the collapse and on the difference 
in velocity Av c between particles in different directions, e.g. 
along the largest and the smallest semi-axis. Given that the 
structure is close to virial equilibrium around t « r c we have 
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Figure 4. Fraction of ejected particles for different values of N. 
The best fit pf = 0.1 + 0.017 log(AT) is a “blind” numerical result, 
while the fit Pf = 0.29 — 2.88 /y/~N will be discussed in Sect |3.21 


Figure 5. Behavior of £8 o(5t c ) and of the initial flattening ratio 
lso (0) (i.e., Eq llH) as a function of N 


2K b « \W b \ and 

A T, c « Uc «vW~^^- (13) 

We thus get 

dt > To) ~ • N 1 ' 6 ~ N - 1/3 (14) 

where we used Eq[l0] Eq[9] and Eafl3l 

We have estimated t for the 80 % more bounded parti¬ 
cles (i.e., t g0 % — this choice avoids the signal to be affected 
by large fluctuations), at tmax = 5r c and we considered 
averages over the different realizations. As shown in Fig[5] 
(where, for comparison, we report the behavior of the initial 
flattening ratio tgo%(0)) the agreement between Eg 1 141 and 
the data is very good over two orders of magnitude in N. 

According to the mechanism of energy gain we have 
discussed above, we expect: the largest semi-axis of (i) the 
initial cloud at t = 0 (ii) the bound particles at tmax and 
(iii) the free particles at t max should be parallel. We call tp 
the angle between (i) and (iii) and ip the angle between (ii) 
and (iiifl. The results are shown in Fig[S] the probability 
density function for both bp and tp are peaked at 0 showing 
that these axes are parallel among each other. In particular, 
ejection occurs more likely in the direction of the the initial 
largest semi-axis, while the final virialized structure is less 
correlated with such a direction. 


4 these are measured from the computation of the corresponding 
eigenvalue A 3 for bound a free particles at t 



cos(angle) 

Figure 6. Probability density function (PDF) for the angles tp 
and p (see text). 

3.2 The ellipsoidal case 

When the initial cloud shape is a prolate ellipsoid then 
its initial asymmetry is amplified by the violent relaxation 
mechanism in a way that is equivalent to that of a spherical 
cloud. However, in this situation the amplitude of the devi¬ 
ation from spherical symmetry, instead of being controlled 
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Figure 7. Black (red) dots show the fraction of particles initially 
with r > R c (r < R c ) that are ejected after the collapse (average 
over 20 realizations). 

by Poisson fluctuations (i.e., Eg flTll is tuned by the ratio 
between the largest 03 and the smallest 02 = ai semi-axis 0. 

Indeed, we find that, as for spherical case, ejected par¬ 
ticles are mostly those which originally lie close to the 
boundaries of the system: i.e. particles for which initially 
ro > R c = ai (see FigO In addition, Fig( 8 ] shows the his¬ 
togram of the values of (see Sect l3.1l for definitions) for 
three different sets of simulations with different initial t( 0 ). 
The result is similar to the spherical case (see Fig( 6 j: the 
largest semi-axis of the initial conditions and of the post 
collapse bound and free particles are parallel among each 
other. 

According to the model of ejection and system defor¬ 
mation discussed for the spherical case, the increase of i( 0 ) 
corresponds to a decrease in the number of particles for an 
initially spherical cloud (see Eg 1111 ) so that: (i) the spread 
in arrival time AT increases and the collapse becomes softer 
as particles arrive at the center at different times, (ii) the 
fraction of ejected particles decreases and (iii) the minimum 
gravitational radius reached during time evolution becomes 
larger: 

• Fig03 shows that the spread in the arrival time AT in¬ 
creases linearly with iso(O). This can be easily explained by 
considering that the spread of arrival time AT is equal to the 
difference between the free fall time respectively of a parti¬ 
cle initially placed at R c ( 1 + i(0)) and of a particle initially 
placed at R c (for small enough i( 0 )): 

AT a {Mr,, (F +-')» f«<») ■ <«> 


Figure 9. Average over 20 realizations of the spread in the arrival 
time AT versus the initial flattening ratio *,(0) for the case of 
ellipsoidal initial conditions (black dots). In addition we plot AT 
versus the actual initial flattening ratio for the simulations of a 
purely spherical cloud with different N (red dots), i.e. *.(0) ~ 
AT-i /2 

— see Fig[5] 

(For the spherical cloud we used Ea llOl and Ea llll to compute 
AT in function of t(0).) 

• Fie llOl shows the fraction of ejected particles as a func¬ 
tion of t 8 o( 0 ) for both the ellipsoidal and the spherical 
clouds: pf decreases with i( 0 ) as the collapse gets softer. 
In particular, we note that pf oc t(0) for both the ellipsoidal 
case at fixed N and varying i(0) > TV -1 / 2 and for the spher¬ 
ical case with varying N and thus with i(0) ~ iV -1 / 2 . This 
result suggests a straightforward interpretation of the frac¬ 
tion of mass ejected in function of the number of particles 
N for the spherical case (see FigUJ) in terms of finite-size 
JV-dependent deformation of the initial spherical cloud. In 
addition, if we extrapolate the behavior for t(0) —7 0, which 
for the spherical case corresponds to N —» oo, we find that 
there is a saturation in the fraction of mass ejected 0 . 

• Fie llll reports the minimum value of the gravitational 

radius during time evolution as a function of t( 0 ) for 

both spherical and ellipsoidal initial conditions, with differ¬ 
ent N and t(0) respectively, together with the theoretical 
fit: 

Kin ~ 40 ) 2/3 (16) 

where we used, for the spherical case, Eq.[9]and Eg II 1 1 Again 
the agreement is extremely good showing that the collapse 
depth is controlled by t( 0 ), which plays, for small enough 
values, the same role of density fluctuations. 

The asymptotic value of tso% as a function of AT x 


5 We consider a situation in which the initial flatness ratio is 
larger than Poisson fluctuations, i.e. i(0) 7§> 1 / 2 


6 For a discuss i on abo ut the Vlasov-Poisson limit of this system 
see I Joyce et alj i2009bh 
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Figure 8. As Fig|6]but for simulations with initial flattening ratio *.(0) = 0.05,0.15,0.25 (from left to right). 



Figure 10. Average over 20 realizations of the fraction of ejected 
particles vs *> 80 (0) (black dots). The case of a uniform sphere with 
different N is also shown (red dots). The dashed line corresponds 
to the fit Pf ~ 1/y/N + const, in FigH] 

R ™ zn , for different initial flattening ratios and averaged over 
20 realizations, is shown in Fig. 1121 We note that the linear 
behavior described by Ea ll2l nicely fit the observed behav¬ 
ior for the spherical simulations. However for the ellipsoidal 
case, when i(0) > 0.1, the asymptotic value of iso clearly 
deviates from the simple growing behavior. This occurs be¬ 
cause of the formation of substructures, a situation that 
clearly becomes more complex than the one we considered 
above. 

In summary the violent relaxation mechanism as illus¬ 
trated by Figs l9llH applies equally well to both the spherical 



Figure 11. Average over 20 realizations of the minimum of the 
gravitational radius vs i(0) (black dots). The case of a uniform 
sphere with different N is also shown (red dots). The dashed line 
represents the fit with Ea ll6l 

and the ellipsoidal initial conditions, considering that finite 
N density fluctuations in the former case play the same role 
of the ratio between the largest and the smallest semi-axis 
in the latter case. 


4 VIRIALIZED AND EJECTED PARTICLES 

We first discuss the properties of the virialized structure 
formed after the gravitational collapse, highlighting that 
some of their main statistical properties, such as the density 
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Figure 12. Evolution of the asymptotic value of l as a func¬ 
tion of AT x aJR™ in , for both ellipsoidal and spherical initial 
conditions. Eq. CCD 


Figure 13. Upper panel: Density profile for the virialized state 
formed in the collapse of a spherical cloud (a line with slope r ~ 4 
is plotted for reference). Bottom panel: particle energy (for bound 
particles only and in normalized units) as a function of the dis¬ 
tance from the center of mass. 


and velocities profiles, are formed during the violent relax¬ 
ation mechanism and do not show an imprint of the initial 
conditions. This result, of course, does not hold in general 
for any cold initial density profile: the identification of the 
phase space parameters of a cold initial cloud for which the 
same density and velocity profiles are formed goes beyond 
the scope of the present work. Then, we describe the dif¬ 
ferent configurations that ejected particles can form, which 
instead depend on the initial cloud shape. 


4.1 The universal density and velocity profiles of 
bound particles 


As mentioned in the introduction, it was recently shown that 
the post-collapse virialized structure, formed when a cloud 
is initially cold enough, has a universal density profile. That 
is, the same profile was found when the initial density profile 
was uniform, power law wit h 0 ^ a ^ 2, Plummer, Gaussian 
or Hernquist JSvlos Labinfeouj ) . The best fit density profile 
is 


P(r) 


Po 

1 + (r/r o) 4 


(17) 


where, in the uniform case, po and ro depend on the num¬ 
ber of points N. This profile was shown to form when bound 
particles can have energies close to zero, so that they are al¬ 
lowed to orbit mostly in radial trajectories around the dense 
core, forming thus the p(r ) ~ r -4 tail (see Fig |13l) . We stress 
that this situation occurs when a fraction of the system par¬ 
ticles get positive energy after the collapse. Thus the mass 
and energy ejection and the formation the p(r) ~ r -4 pro¬ 
file are two aspects of the same phenomenon intrinsic to the 
violent relaxation mechanism. 


Fig |14l shows that Eg 1 171 gives an excellent fit also when 
the initial condition is an ellipsoid; this fact suggests a com¬ 
mon origin for the density profile as an outcome of the vi¬ 
olent relaxation mechanism. In addition, the radial veloc¬ 
ity dispersion shows a Keplerian behavior for all cases (see 
Fig l 161) which is a further key element in the physical model 
for the formation of the ~ r~ 4 tail. This means that the ini¬ 
tially ellipsoidal cloud, in the range of parameter space con¬ 
sidered here, shows the same gravitational collapse dynamics 
as the spherical clouds. In particular, the virialized structure 
emerging from the violent relaxati on share the s ame density 
and velocity profiles. As argued in ISvlos Labinil (120121 . 120131 ) 
these are formed during the violent relaxation mechanism 
because of the large energy changes which particles undergo 
during the collapse. 

As noticed above, the remnants of cold collapses of ini¬ 
tially ellipsoidal clouds are strongly non-spherical. In order 
to take into account the elongated cloud shape, ellipsoidal 
shells of particles are used in the density profile calculation: 
in particular, we computed the density profile in ellipsoidal 
shells with parameters proportional to (and smaller of) those 
of the final structure. In Fig llSI we show density profiles and 
their respective fits, using elliptical and spherical shells for 
the case i(0) = 0.15. While the free parameters po and ro 
are different for these two profiles, one can see that Eq |17l 
still provides an excellent fit to the density profile. 

The distinctive feature of the violent relaxation mech¬ 
anism is represented by the large particle energy change. 
This is illustrated in Fig llTl (lower panel), where it is plot¬ 
ted the initial and final particle energy distribution P(e) for 
t(0) = 0.05 and i(0) = 0.25. As the collapse is softer in the 
former case, the final P(e) has a smaller spread than for case 
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Figure 14. Upper Panel: Density profile for clouds with differ¬ 
ent initial shape: spherical and prolate ellipsoidal with i(0) as 
reported in the labels. As a reference we have also reported 
the behavior of the quasi equilibrium profile. Lower panel: par¬ 
ticle distribution energy at both the initial and the final time 
for ellipsoidal initial conditions respectively with i(0) = 0.05 and 
i(0) = 0.25. The particle energy e has been normalized to eo, the 
absolute value of the the potential energy of a particle placed at 
distance R c from the center of a uniform spherical cloud. 



Figure 15. Profiles calculated with ellipsoidal shells and spheri¬ 
cal shells, for the case i(0) = 0.15. 


4.2 The asymmetric feature of particle ejection 

As discussed above, because ejected particles are those ini¬ 
tially at ai = d 2 ^ r St < 23 , they are not spherically symmet¬ 
ric neither initially nor for t > r c : indeed, the collapse ampli¬ 
fies their initial asymmetry resulting in a very anisotropic fi¬ 
nal distribution. To characterize their shapes we measure the 
t, r and (f> (see the upper panel of Fig l 171) : we may see that by 
increasing t( 0 ), ejected particles do not maintain spherical 
symmetry (as for i(0) = 0). Instead, they first tend to form 
a flat structure and then a bar for i(0) >0.1 (see Fia ilTI) . 
Obviously particles ejected on the two opposite sides of the 
virialized structure will have velocities oriented in opposite 
directions, i.e. they will be anti-correlated. When the initial 
ellipsoid is oblate then ejected particles form a triaxial struc¬ 
ture. Indeed, Fig llSI (upper panel) shows that 1 « 2 <j> and 
t « 1/2. Finally for a triaxial initial cloud ejected particles 
form an anisotropic structure that clearly depends on the 
choice of the initial axes (see Fig |19l) . 

As we can see in Figs. 11711191 and [IS] by changing the 
value of the initial flatness ratio i( 0 ) and the initial shape of 
the structure (i.e. prolate, oblate or triaxial) one can form 
different final distributions of ejected particles which are 
close to flat. While in most cases the majority of the free 
particles are ejected within a bar along the major semi-axis 
of the final structure, we note that for the oblate case, espe¬ 
cially when i(0) = 0.05,0.1, the distribution of the ejected 
particles is closer to that of a disk. 


t(0) = 0.05: in particular, both the small and high energy 
tails are reduced. This corresponds to the fact that the final 
core is less dense and there are fewer ejected particles. 


5 DISCUSSION AND CONCLUSIONS 

The gravitational collapse of an initially cold and ellipsoidal 
cloud of particles shows non-trivial characteristics: some of 
them are similar to the ones observed when an initially 
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Figure 16. Radial velocity dispersion for clouds with different 
initial shape. As a reference we have reported also the Keplerian 
1 jr behavior. 


spherical cloud is considered as an initial condition. In par¬ 
ticular, we have focused our attention on the the formation 
of power-law density and radial velocity profiles and on the 
ejection of a fraction of mass and energy from the system. 
The analogy between the spherical and the ellipsoidal cloud 
collapse can be understood considering that the initial de¬ 
viations from spherical symmetry in the former case plays 
the same dynamical role of density fluctuations in the latter 
one. 

Indeed, we have shown that in the ellipsoidal case the 
ratio between the largest as and smallest ai semi-axis play 
the same role of density fluctuations for a spherical cloud. 
Indeed, N randomly placed particles in a sphere, have an in¬ 
trinsic ellipsoidal deformation oc y/N due to Poisson fluc¬ 
tuations; thus increasing N has the same effect of decreasing 
the difference between as and a i. We have discussed in de¬ 
tail that the mechanism of particles energy gain during the 
collapse gives rise to the ejection of those particles that are 
initially placed close to the cloud boundaries: for an ellip¬ 
soidal cloud this implies an asymmetric ejection. 

Thus, the virialized state formed after the collapse of 
an ellipsoidal cloud shows the same features of the cold 
spherical cloud remnant: in both cases the density and ra¬ 
dial velocity profiles arise from the violent relaxation mech¬ 
anism. Spherical clouds with different density profiles can 
however form different virialized structures and a more com¬ 
plete study of the cloud parameters phase space is necessary 
to draw more general conclusions about quasi-equilibrium 
profiles. 

Note that the mechanism described here is, in principle, 
dif ferent from t he the radial orbit instab il ity discussed by 
lAntonovI (ll96ll ): lFridman fc Polvachenkol (1 1984 ). This insta¬ 
bility characterizes a different case from the one considered 


3 

2 
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Figure 17. Upper Panel: i,t, </> of ejected particles for a prolate 
initial ellipsoidal cloud as a function of the initial i(0). Middle 
and lower panels: Projections onto two different planes of both 
ejected and bound particles (that are clustered in the center) for 
prolate initial condition with N = 10* and /,(()) = 0.25. The bar¬ 
like shape of ejected particles is clearly visible. 
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Figure 18. As Fig|l7lbut for oblate initial clouds with N = 10 5 
and t(0) = 0.05. 


Figure 19. As Fig U7l but for triaxial initial clouds with N = 10 5 
with axis as = 1.1, a 2 = 1.05, ai = 1. 
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here, i.e. a spherically symmetric stationary solutions of the 
collisionless Boltzmann equation with purely radial orbits. 
The radial orbits instability characterizing equilibrium mod¬ 
els and systems undergoing to strong collapses and violent 
relaxations are distinct, but probably c onnected, physical 
mecha nisms. Indeed, as firstly noticed bv I Merritt fe Aguilar! 
11198511 a radial instability can be the mechanism that, in 
the collapse of an initially cold system in which particles 
have mainly radial orbits, gives rise to a triaxial struc¬ 
ture. Our findings clarify the dynamics driving the forma¬ 
tion of a bar/disk structure in more detail than in previ¬ 
ous works, by focusing on the relation between initial and 
final conditions for simple non-spherical systems. However 
a detailed theoretical understanding of the radial orbit in¬ 
stability for non-equilibrium collapsing systems is still an 
open problem: we refer the interested reader to a forthcom- 
ing wo rk for a more deta il ed discussi on about this point 
(ISvlos Labini. Benhaiem fe Jovcell2014l 'l. 

On the other hand ejected particles form non-trivial 
anisotropic configurations, like bars or disks, that depend 
on the initial cloud shape: a small initial deviation from 
spherical symmetry is amplified by the violent relaxation en¬ 
ergy gain mechanism and gives rise to highly non-spherical 
clouds. The ellipsoidal clouds that we have considered rep¬ 
resent thus a further step toward a more realistic cloud col¬ 
lapse, that has allowed us to perform controlled N-body ex¬ 
periments to quantify the effect of the initial sphericity. 

This result, although obtained for a very simplified and 
unrealistic physical cloud model, is interesting as it points 
out that it is relatively natural to generate, from an initial 
condition close to isotropic, a central (host) structure sur¬ 
rounded by particles, or even clusters of particles, that are 
strongly anisotropically distributed, even sometimes close 
to a planar distribution. Galaxy formation is a much more 
complex process, involving physical mechanisms such as 
gas dynamics, effects of dark matter halos, etc., which are 
not taken into account in this work. However, it is worth 
mentioning that recently several observational studies have 
pointed out that dwarf satellites of galaxies both the Milky 
Way and the Andromeda galaxy are no n isotropica ll y dis¬ 
tribu ted around their host galaxy (see llbata et al.l (120131 . 
120141 1 and references therein). In particular it was observed 
that dwarf satellites are aligned in a thin plan structure and 
that they have coherent kinematic properties. Our results 
suggest that these particular features might be better un¬ 
derstood from the collapse of a non-spherical cloud. In a 
forthcoming work we will consider more realistic initial con¬ 
ditions and a detailed comparison with observations. 

We thank Michael Joyce for useful discussions and com¬ 
ments. We also thank an anonymous referee for a number of 
interesting comments that have allowed us to improve the 
presentation of our work. Numerical simulations have been 
run on the Cineca supercomputer (project ISCRA QSS- 
SSG) and on the super-computer Mesu of UPMC. 
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